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Introduction 


The  objective  of  this  effort  is  to  develop,  implement  and  evaluate  signal  generation 
and  processing  techniques  oriented  towards  improving  the  detection  and  estimation  of  targets 
and  other  parameters  in  noisy  and  cluttered  environments.  It  is  also  desired  to  quantify  any 
improvements  in  the  detection  and  estimation  of  targets  in  noisy  and  cluttered  environments. 

In  developing  new  or  enhancing  existing  signal  processing  algorithms  for  improved 
detection  and  estimation  performance  of  exiting  and  planned  sensors,  one  has  to  be  aware  of 
environmental  induced  effects  and  how  the  algorithm  optimizes  the  sensor  performance.  In 
many  cases,  optimization  is  achieved  by  generating  a  set  of  weights  and  then  applying  them 
to  the  elements  of  both  array  elements  and  time  delay  elements  of  the  system.  Space  Time 
Adaptive  Processing  (STAP)  is  one  example  of  these  adaptive  processes.  Performance 
evaluation  and  a  fair  comparative  analysis  of  these  techniques  require  the  developer/user  to 
have  a  signal  generation  that  is  compatible  with  real  world  data.  However,  modeling  these 
effects  is  a  very  tedious  and  difficult  task.  For  these  reasons,  alternative  approaches  to  the 
physical  model  signal  generation  are  sought.  This  will  speed  up  analyses.  Furthermore 
Monte  Carlo  type  analyses  can  be  done  in  a  timely  manner.  The  “representative  model  that 
was  suggested  by  Dr.  James  Michels  of  the  Air  Force  Research  Laboratory  (AFRL),  Sensor 

Directorate  (SNRT)  is  one  such  example  [1]. 

In  this  effort,  we  discuss  three  methods  for  signal  generation.  First,  the  representative 
model  of  Dr.  James  Michels  is  introduced.  Next,  the  physical  model  developed  by  Scientific 
Studies  Corporation  (SSC)  [2]  will  be  briefly  discussed.  We  will  refer  to  it  in  the  remainder 
of  this  report  as  the  physical  model.  Finally,  a  third  technique,  which  makes  use  of  both  the 
physical  model  and  the  representative  model,  is  introduced  [3].  This  approach  will  be 
referred  to  as  the  hybrid  model.  This  set  of  data  signal  generation  forms  the  basis  of  our 
simulation  capability  tool. 

This  data  simulation  tool  is  capable  of  generating  time  series  data  as  viewed  by  a 
linear  phased  array  on  an  airborne  moving  platform.  To  avoid  ambiguities  and  for 
completeness  of  the  work,  the  same  notation  as  that  developed  by  SSC  will  be  used  in  this 
report.  Using  the  representative  model  for  example,  the  idea  is  to  generate  data  directly  from 
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a  specified  covariance  matrix  for  evaluations  of  STAP  algorithms.  The  method  provides 
control  of  both  spatial  and  temporal  correlation  as  well  as  the  probability  density  function 
(pdf)  of  the  processes.  Parameters  such  as  partial  temporal  and  spatial  correlation  of  signal 
and  interference  resulting  from  the  observation  platform  motion,  beam  pointing  and  internal 
clutter  motion  are  also  included. 

This  report  is  organized  as  follows.  Chapter  2  describes  the  models  for  signal 
generation  used  in  this  effort.  Examples  involving  the  representative  model  will  be 
introduced  in  section  3.  True  covariance  computation  and  time  series  data  generation  are 
discussed  in  Chapter  4.  In  chapter  5,  several  techniques  for  data  generation  will  be 
discussed.  In  chapter  6  we  show  some  examples  using  the  proposed  models.  Concluding 
remarks  and  recommendation  for  future  work  will  be  discussed  in  chapter  7. 
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2 

Phased  Array  Radar  Model 


In  this  section,  a  system  model  is  developed.  This  system  is  commonly  used  for 
airborne  surveillance  pulsed-Doppler  radars.  The  model  is  based  on  a  side  looking  phased 
array  radar  configuration  with  a  scanning  capability  of  ±90°.  Note  that  each  component 
(target,  jammer,  clutter  and  noise)  is  modeled  independently  of  the  others.  It  is  our  goal  and 
intent  to  generate  true  covariance  matrices  for  each  of  these  components  as  well  as  to 
generate  a  sequence  of  space-time  independent  data. 

Targets  are  modeled  as  point  targets  in  azimuth  and  elevation  and  for  a  wide  range  on 
scenarios.  This  includes  the  constant  amplitude  case  and  Swerling  cases  1  through  4. 
Jamming  signals  are  modeled  as  barrage  type  noise.  This  assumes  point  sources  with 
Gaussian  distributed  wide-band  noise  in  both  azimuth  and  elevation.  Receiver  noise  is 
assumed  to  be  uncorrelated  in  both  time  and  space.  Ground  clutter  is  modeled  as  a  large 
number  of  statistically  independent  clutter  patches  with  each  patch  consisting  of  a  large 
number  of  individual  re-radiators. 

The  array  is  assumed  to  be  linear  uniformly  spaced  with  element  spacing  denoted  by 
d.  It  consists  of  J  identical  elements  (or  identical  beamformed  sub-arrays).  The  platform  is 
assumed  to  be  moving  at  a  velocity  vp  at  an  altitude  Hp.  The  array  geometry  is  shown  in 
Figure  1. 
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Figure  1.  Airborne  Surveillance  Scenario 


The  antenna  boresight  is  assumed  to  be  normal  to  the  platform  direction.  However, 
the  platform  velocity  vector  vp  may  have  an  offset  angle  from  the  platform  direction.  This 
results  in  a  crab  angle  denoted  by  5. 
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Figure  2.  Airborne  surveillance  geometry 
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Representative  Model 


The  approach  consists  of  first  forming  a  two-dimensional  (2-D)  spatial  and  temporal 
correlation  function  with  a  corresponding  2-D  spectrum  as  shown  below. 


Figure  3.  Correlation  function 
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The  term  “ representative ”  is  used  to  denote  the  generic  shaping  of  this  correlation 
and  its  corresponding  spectrum  through  the  use  of  parametric  functions  referred  to  as 
“ shaping ”  functions.  This  procedure,  unlike  the  physical  model,  does  not  account  for  a 
detailed  physical  description  that  would  then  give  rise  to  such  a  function.  This  approach 
shapes  the  correlation  function  with  the  intent  to  obtain  profiles  representative  of  those 
observed  in  practice.  Figure  5  and  6  show  the  correlation  function  and  its  corresponding 
spectrum,  respectively,  for  the  case  when  the  radar  platform  is  stationary  and  the  mainbeam 
is  pointed  at  0°  (boresight). 
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Figure  5.  Correlation  function,  stationary  platform,  <j>  =  0°. 


Power  Spectrum 


Doppler  Frequency  ~°  5  ~0  5 _ Sgatiamgguenc^ 

Figure  6.  Corresponding  spectrum. 


Note  from  the  above  figures  that  the  Doppler  spread  is  due  to  internal  clutter  motion. 
The  spatial  distribution,  on  the  other  hand,  is  the  result  of  the  antenna  beam  profile. 


To  show  the  simplicity  of  the  technique,  assume  that  the  platform  is  still  stationary, 
but  the  mainbeam  is  pointed  at  and  angle  (j)  =  20°.  This  results  in  the  following  spectrum  in 
the  angle-Doppler  domain. 


Figure  7.  Spectrum,  vp  =  0,  $  =  20°. 


Assume  now  that  the  platform  velocity  is  non-zero.  This  result  in  a  Doppler-spatial 
frequency  dependence  that  yields  the  rotated  spectrum  (see  Figure  8  below). 


The  shaping  functions  described  above  are  used  to  control  both  the  temporal  and  the 
spatial  profiles  of  the  correlation  distributions.  The  values  of  the  cross-  and  auto-correlation 
functions  are  then  used  as  the  elements  of  the  multi-channel  correlation  matrix.  This  matrix 
is  then  used  in  the  process  synthesis  procedures  to  generate  time  series  data  that  contain  the 
correlation  specified  in  the  matrix. 
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4 

Correlation  Matrix  Description 


Define  the  zero-mean  (Jxl)  observation  data  vector  X(n)  at  time  n,  under  the  null 
and  alternative  hypothesis  H0  and  H, ,  respectively,  as 

H0 :  X(n)=C(n)+l(n)+W(n)  ;  n=l,2,---,N  (1) 

H, :  X(n)=S(n)+C(n)+l(n)+W(n)  ;  n=l,2,---,N  (2) 

where  S(n)  is  the  signal  vector,  C(n)  is  the  temporally  correlated  clutter  process,  l(n)  is  the 
interference  process  and  W(n)  is  the  additive  white  noise  process  assumed  to  be 
uncorrelated  in  both  space  and  time. 

For  a  fully  adaptive  signal  processor,  a  (JNxl)  concatenated  observation  data  vector 
X  is  formed  as 

w 

x=  x'2) .  o: 

,X(N). 

4.1  Signal  Model 

In  this  section,  we  review  the  different  signal  models  employed  in  our  simulations. 


4.1.1  Constant  Amplitude  Signal 

The  signal  model  provides  a  (Jxl)  vector  time  sequence  S(n).  For  the  constant 
amplitude  signal  model,  we  have 

S(n)=aU(n)  (4) 

where  a  represents  the  complex  signal  amplitude  and  the  signal  steering  vector  U(n)  is  given 
by 

U(n)=Uexp{j2n(n-l)ti3t}.  (5) 

Note  from  the  above  Equation  that  the  spatial  steering  vector  U  is  defined  to  be 


11 


u= 


(6) 


1 

exp{j27t$,} 

exp{j27c(j-l)dt} 


In  Equation  (5)  and  (6),  G5t  and  denote  the  normalized  target  Doppler  and  spatial 
frequencies  as  viewed  from  the  platform,  respectively,  and  are  given  by 


G3t  =T{ft  +fp  cos(et  )sin(4>t  -  y)} 

(7) 

and 

tit=(d/X)cos(et)sin(<|>t  -y). 

(8) 

In  the  above  Equations,  the  frequencies  ft  and  fp  are  given  by 

f,=(2v,A.) 

(9) 

and 

f,=(2v,AJ. 

(10) 

The  complex  signal  amplitude  a  can  be  expressed  as 

a=Aexp(j80) 

(ID 

where  A  is  the  magnitude  and  80  is  the  initial  phase.  Therefore  the  k-th  channel  signal  is 
expressed  as 

s  k  (n)=  A  exp(j80  )exp{j2rc(k  -  l)dt  }exp{j27t(n  -  l)tiJ  t } ;  n = 1, 2,  •  •  • ,  N ,  k = 1, 2,  •  •  • ,  J .  (12) 


4.1.2 

Case  1:  Constant  Amplitude,  Random  Phase  Signal  Model 

In  the  constant  amplitude,  random  phase  signal  model,  A  is  kept  constant  over  both 
the  time  pulses  of  the  CPU  as  well  as  over  multiple  ensemble  realizations  (i.e.,  over  multiple 
CPI’s).  The  phase  80  is  kept  constant  over  the  time  pulses  of  the  CPI,  but  is  assumed  to  be  a 

random  variable  uniformly  distributed  in  the  interval  (0,27t),  over  multiple  ensemble 
realizations. 
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4.1.3 

Case  2:  Constant  Amplitude,  Constant  Phase  Signal  Model 

In  the  constant  amplitude,  constant  phase,  signal  model,  both  A  and  80  are  kept 
constant  over  both  the  time  pulses  of  the  CPU  as  well  as  over  multiple  ensemble  realizations 
(i.e.,  over  multiple  CPFs. 

4.1.4 

Case  3:  Swerling  I  Model 

The  Swerling  I  model  characterizes  a  signal  with  perfect  correlation  pulse-to-pulse 
over  a  coherent  processing  interval  (CPI)  consisting  of  N  pulses,  but  random  over  successive 
realizations.  In  this  case,  the  amplitude  A  is  modeled  as  a  Rayleigh  distributed  random 
variable  across  realizations  but  kept  constant  over  a  given  CPI.  Its  mean  value  is  A.  The 
phase  80  is  also  kept  constant  over  a  single  realization  of  N  time  samples  (one  CPI),  but  is  a 

random  variable  uniformly  distributed  between  (0,2rc),  over  multiple  ensemble  realizations. 

4.1.5 

Case  4:  Swerling  II  Model 

The  Swerling  II  model  characterizes  a  signal  that  is  totally  uncorrelated  from  pulse  to 
pulse  over  a  coherent  processing  interval  (CPI)  consisting  of  N  pulses.  In  addition,  it  is  also 
assumed  to  be  a  random  variable  over  successive  realizations.  The  amplitude  A  is  modeled 
as  a  Rayleigh  distributed  random  variable  across  realizations  but  kept  constant  over  a  given 
CPI.  Its  mean  value  is  A.  The  phase  80  is  also  assumed  to  be  a  random  variable  uniformly 

distributed  between  (0,2tt),  over  multiple  ensemble  realizations. 

5.1.6 

Case  5:  Random  Signal  Model,  Partial  Temporal  Correlation 

The  random  signal  model  contains  temporal  correlation  over  the  N  time  samples  and 
is  uncorrelated  over  multiple  realizations. 

4.2  Interference  Model 

The  partially  correlated  noise  jamming  signal  model  is  defined  as 

^  (n ,  0J  )=  l(n)  •  exp{j2  Jt(m  -  l)fy  }•  exp[j7t(n  - 1)05)  ]  (13) 
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where  l(n)  is  a  complex  noise  process  with  partial  temporal  correlation,  d;  is  the  normalized 
interference  source  spatial  frequency  and  05;  is  the  normalized  interference  center  Doppler 
frequency.  The  quantities  and  05;  can  be  expressed  as 

■d;  =(d/A,)cos(0;  )sin(<J>; )  (14) 

and 

05;  =(vi/X)r+(2vp/X.)rcos(0i)sin(<|)i),  (15) 

where  vp  is  the  platform  velocity  and  V;  is  the  interference  source  velocity.  For  white  noise 

barrage  jamming  signals,  we  have 

^(n,ei)=l(n)*exp{j27t(m-l)di}  (16) 

where  l(n)  is  a  complex  white  noise  process. 


4.3  Clutter  Model 

In  the  following,  three  (3)  procedures  for  the  determination  of  the  known  clutter 
covariance  matrix  applicable  for  the  airborne  radar  problem  are  described.  Specifically,  the 
auto-  and  cross-correlation  functions  for  the  clutter,  signal  and  interference  sources  as  viewed 
from  airborne  phased  array  radars  are  developed.  These  functions  will  then  be  used  to  form 
the  known  covariance  matrices  for  data  synthesis.  These  include  (1)  the  representative 
model,  (2)  the  physical  model  and  (3)  the  hybrid  model.  These  methods  lead  to  a  procedure 
for  the  generation  of  a  clutter  time  series  with  correlation  provided  by  the  known  covariance 
matrix.  Note  that  Dr.  James  Michels  has  suggested  both  the  representative  model  and  the 
hybrid  model.  Scientific  Studies  Corporation  Inc.  (SSC)  has  coded  the  physical  model. 
Minor  changes  have  been  made  to  this  code  to  fit  both  the  Rome  Laboratory  Space-Time 
Adaptive  Processing  (RLSTAP)  needs  and  our  own  needs. 

In  this  section,  the  expressions  for  the  auto-  and  cross-correlation  functions  between 
side  mounted  airborne  sub-arrays  are  derived.  Note  that  only  azimuth  variations  within  a 
given  range  resolution  cell  are  considered.  Two  approaches  are  discussed.  The  first 
approach  does  not  take  into  account  the  effect  of  platform  crab  angle  and  the  second 
considers  this  effect.  First,  we  do  not  consider  the  effects  of  platform  crab  angle. 
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Let  ck  (n)  denote  the  complex  baseband  clutter  process  from  the  k-th  channel  at  time 
n.  Assuming  stationarity,  the  clutter  cross-  channel  correlation  function  between  channels  k 
and  i  is  given  by 

Rc,ki  (l,  )=Efe  (n)c*  (n  -  lt  )}=E{ck  (n  +  lt  )c*  (n)}, 
where  1,  is  the  temporal  lag  value.  A  functional  form  for  the  correlation  function  Rc>ki(lt)  is 
now  proposed.  To  this  end,  let  <|>0  be  the  azimuth  angle  to  the  center  of  the  beam  position  as 
measured  from  boresight  of  the  antenna.  The  cross-correlation  function  between  the  k-th  and 
i-th  channels  is  expressed  as 

Rckj  (lt  )=ac,kac,i  Fc  (lt,k,i)exp|j27t(k  -i)|sin(<))0)|exp|j2rftT 

(-N  +  1)  <lt  <(N-1) 

where 

•  ac  k  and  ac .  are  the  standard  deviations  of  the  clutter  process  as  observed  on 

channels  k  and  i,  respectively 

•  lt  is  the  temporal  lag 

•  Fc  (lt  ,ls )  is  a  two-dimensional  function  of  lt  and  ls  for  the  clutter  which  shapes  the 
magnitude  of  the  clutter  cross-correlation  function  in  both  the  temporal  and  spatial 
domain,  also  referred  to  as  the  shaping  function 

•  T  is  the  inter-pulse  period 

•  vp  is  the  platform  velocity 

•  X=  (c/f o )  is  the  transmitter  wavelength,  f0  is  the  transmitter  frequency  and  c  is  the 
speed  of  light 

•  d  is  the  inter-element  spacing 

•  (|)0  is  the  center  of  the  beam  position  in  azimuth. 

Let  ls  be  the  spatial  lag,  i.e.;  ls  =(k  -i),  then  Equation  (17)  becomes 

R  e  (1 1 ,  ls  )= aCikaCiiFc  (1, ,  ls  )exp|j27tls  ^  sin(4>0  )j  exp  j  j2idtT  -^sinfo0  )|  ( j 

(-N  +  l)<lt£(N-l)  ;  (-M  + 1)  <ls  < (M - 1) 
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The  two-dimensional  (2-D)  power  spectral  density  S(f,,fs)  is  obtained  by  taking  the  2-D 
FFT  of  Rc  (l,  ,ls ).  The  frequencies  ft  and  fs  are  referred  to  as  the  temporal  (Doppler)  and 
spatial  frequencies  and  are  given  by 

f,=^sin(<l>0)  (19) 

and 

f^sinfao),  (20) 


respectively. 

Using  Equations  (19)  and  (20),  it  can  be  easily  shown  that 


Let  tan  (a)  be  the  quantity 


tan(a)= 


(21) 


(22) 


Then  ft  and  fs  are  related  through  the  equation  ft=tan(a)f$.  Note  that  tan  (a)  is  the  slope 
of  the  clutter  in  the  normalized  Doppler-spatial  frequency  domain. 

Therefore  a  is  given  by 


a=tan 


(23) 


\  J 

For  a  non-zero  platform  crab  angle,  yc ,  sin(<])0 )  in  the  above  equations  is  replaced  by 
sin(<|)0  +  \|/c).  This  results  in  a  non-linear  relation  between  ft  and  fs  in  Equation  (23).  This 
results  in  a  power  spectrum  that  is  no  longer  centered  along  the  linear  slope  in  the  normalized 
Doppler-spatial  frequency  domain,  but  rather,  along  the  arced  curve.  This  case  will 
discussed  further  in  the  next  section. 

Next,  we  introduce  the  coordinates  (lA,lB),  in  which  the  correlation  function  remains 
separable  in  time  and  space.  The  2-D  shaping  function  is  expressed  as 


(24) 


where  the  corresponding  rotation  in  the  correlation  space  is  achieved  using 
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1A  _  cos(a)  -sin(a)Yl, 

1B  sin  (a)  cos(a)J|_ls_ 

Then,  Equation  (18)  is  rewritten  as 

Fc(lt,ls)=Fct{lt  cos(a)-ls  sin(a)]Fcs{lI  sin(a)+ls  cos(a)}. 

If  the  platform  is  not  moving,  vp  =0 ,  therefore  oc=0 ,  Equation  (26)  becomes 

Fc(lt,ls)=Fc>t(lt)FCiS(ls). 

This  shows  that  Fc(lt,ls)  is  separable. 

4.3.1  Temporal  Shaping  Functions 

Temporal  shaping  functions  can  take  on  several  distributions.  For  a  Gaussian  shaped 
clutter  power  spectrum,  the  temporal  function  Fct(»)  is  also  Gaussian.  In  this  case,  Fc  t(lA) 

can  then  be  expressed  as 

<28> 

The  parameter  pc  t  is  the  one-lag  temporal  correlation  clutter  parameter  that  controls  the 
width  of  the  correlation  function  magnitude  along  the  1A  axis  and  has  range  0<  jxc  t  <  1 .  The 
parameter  1A  mi  is  the  temporal  lag  at  which  the  magnitude  of  the  clutter  cross-correlation 

function  achieves  its  peak  value.  For  the  phased  array  radar  problem  considered  here,  the 
parameter  1A  mi  is  negligible.  It  is  important  for  cases  where  receiver  elements  are  widely 

spaced.  In  this  case  the  time  delay  between  channels  is  considerable.  Only  the  case  of 
1A  ^  =0  is  considered  in  this  section. 

In  the  case  of  a  Gaussian  shaped  function,  the  parameter  |ic  t  is  expressed  as 

m,=exp{-27tS!T!},  (29) 

where  E,2  is  the  variance  of  the  Gaussian  shaped  power  spectrum  measured  along  the 
temporal  frequency  axis  ft .  Note  that  the  form  of  Equation  (18)  will  change  for  other 
shaping  functions.  In  the  case  of  an  exponentially  shaped  function,  Fct(*)  is  expressed  as 

Fc,t(lA)=k,)'Al  (30) 


(25) 

(26) 

(27) 
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where 


Therefore,  Fc ,  (lA )  is  given  by 


|Ac.  =exp{—  2^T }. 


F  (,  )=_«! _ 

AC,t\1A/  2  (n,  1  \2  * 

a  +(2tc1a) 


(31) 


(32) 


4.3.2  Spatial  Shaping  Functions 

The  spatial  shaping  functions  for  the  clutter  Fc  s(ls )  is  expressed  as 

0.  )=lGT  (1,  T 1° .  (<t>0  >3:  (<)>„  ]| ,  (33) 

which  can  also  be  written  as 

Ft'(l,HOT(UW(*„)-  (34) 

In  Equations  (33)  and  (34),  the  term  |Gt  (ls  )|2  is  the  transmitter  magnitude  gain  term  in  the 
correlation  domain  whereas  the  complex  receive  gain  terms  Gk  (<])0)  and  G;  (<|>0)  for  channels 
k  and  i,  respectively,  are  assumed  to  have  uniform  unit  amplitude  gain  on  each  channel.  The 
term  cos(<j)0 )  represents  the  individual  element  receiver  gain  for  a  dipole  pattern.  The  term 

|Gt  (ls  f  can  be  expressed  as 


|GX  (ls  f  =|p  +  (i  - 


(-M  +  l)<ls<(M-l). 


(35) 


Note  in  the  above  Equation  that  it  is  possible  to  implement  the  Hanning  ((3=0.5),  Hamming 
((3=0.543478261)  and  Blackman-Harris  ((3=0.53856)  element  weighted  gain  profiles. 

In  the  (l  A ,  1 B )  coordinate  system,  we  have 

,  11 2 

(36) 


J  [3  i  ( 1 

L  *.  .  11 

<  p  +  \i  —  pjcos 

(N  -  l)sin(a)+  (M  -  l)cos(a)JJ 

for 

-  [(N  -  l)sin(a)+  (M  -  l)cos(a)]<  1B  <  [(N  -  l)sin(a)+  (M  -  l)cos(a)]. 
For  a  Gaussian  shaped  spatial  function,  we  have 
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Fc.s  Ob  )=  (^c.s  )°B )2  cos2  (<l>o  ) 

where  (ji.cs )  is  defined  as  the  one-lag  correlation  clutter  parameter  and  has  range  (o<|ics  <l). 
The  parameter  (jics )  has  the  same  role  as  (jic  t ).  Whereas  (jic  t )  controls  the  width  of  the 
temporal  frequency  spectrum,  (jic  s )  controls  the  width  of  the  spatial  frequency  spectrum 
related  to  the  beam  width. 

4.3.3  Correlation  Functions 

The  cross-correlation  function  between  channels  k  and  i  is  given  by 
RCJd(lt,ls)=oc,kaciFct{lt  cos(oc)-ls  sin(<x)]FCiS{lt  sin(a)+ls  cos(a)} 

f  2vp 

expj  j27iltT  — 

The  first  complex  exponential  contains  the  phase  delay  between  the  elements  from  a  source 
with  phase  coming  from  <|) =<J>0 .  This  results  in  a  peak  in  the  spatial  spectrum  at 

f,=T|sin(*0).  (38) 

The  second  complex  exponential  term  contains  the  pulse  to  pulse  phase  delay  due  to  the 
platform  velocity.  This  again  results  in  a  peak  at  ft  =fdc ,  where  fdc  is  given  by 

fdc=T^sin((t>0).  (39) 

Note  that  f  dc  is  the  center  of  the  Doppler  shifted  clutter  spectrum. 

The  terms  ac  k  and  oc  i  control  the  magnitude  of  the  cross-correlation  function  that 
itself  controlled  by  the  temporal  and  spatial  functions  Ft(»)  and  Fs(»),  respectively.  As  was 
mentioned  before,  the  width  of  the  spectrum,  along  the  temporal  frequency  f  t ,  is  due  to  the 
internal  motion  of  the  clutter,  described  by  Ft  (•),  and  the  platform  motion  induced  Doppler 
spread  from  all  azimuth  angles  viewed  by  the  antenna  beam  pattern  and  is  described  by 
Fs  (•).  The  latter  is  the  dominant  contributor  to  the  width  of  the  spectrum  in  the 
angle/Doppler  domain.  Note  that  a  narrow  width  in  the  correlation  domain  corresponds  to  a 
wide  spectral  width.  Since  the  width  in  the  correlation  function  Rc  (l,  ,ls )  along  the  1,  and  ls 


•  expj  j27tls  -sin(<|>0)J 
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axes,  is  determined  by  the  product  F,  (•)  Fs  (•),  the  narrowest  of  these  two  functions  in  the 

correlation  domain  yields  the  large  spectral  width. 

Consider  now  the  auto-correlation  function;  i.e.,  k  =  i.  In  this  case,  Equation  (37) 


becomes 

R,Ut(lnls)=<kFc.,{l.  cos(a)-l,  sin(a)}Fcs{lt  sin(a)+ls  cos(a)} 
•  exp|j27tls  ^sin(4>0  )|exp|j27tl,T  -^-sin^  )| 

To  further  simplify  this  equation,  let  ls  =k  —  i=0.  In  this  case, 

exp|j2jds^sin(<t)o)|=l, 


and  Equation  (40)  becomes 


2v„ 


(40) 


(41) 


RCikk0t  )=cJc,kFc,t{lt  cos(a)]Fcs{lt  sin(ot)}expj  j2rtltT  — Lsin(([>0)f .  (42) 


Note  from  Equation  (42)  that  the  auto-correlation  function  has  no  dependence  on  ls.  Assume 
now  that  vp  =0 ,  this  implies  that  a=0 .  The  function  Fe  s{l,  sin(a)}  is  unity  and  the  function 
Fct{lt  cos(a)}  reduces  to  Fc ,  (lt ).  In  this  case,  the  width  of  the  auto-correlation  function  is 
only  dependent  on  the  internal  clutter  motion.  For  the  case  where  1,  =0 ,  Equation  (42) 
reduces  to 

R,*(i.)=R«j,0i)=<..  <43> 

which  is  the  variance  of  the  clutter  process.  Equations  (42)  and  (43)  provide  analytic 
expressions  for  the  elements  of  the  clutter  covariance  matrix  for  the  k-th  range  Rc  k .  In 

section  5,  this  matrix  will  be  used  to  develop  several  procedures  for  data  synthesis. 


4.3.4  Physical  Model  Clutter  Covariance  Matrix 

In  this  section,  a  brief  description  of  the  SSC  physical  model  will  be  given.  For  a 
complete  description  of  the  technique,  please  see  [].  The  equations  that  will  be  described 
here  are  provided  here  for  convenience  only. 

Let  the  JxJ  clutter  covariance  matrix  be  given  by 
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(44) 


N  -1 

R c  (m)=  E|c(n)bH  (n  - m)]=  £ <  exp{j2ranG5p  ^  (dp  )b”  (dp ) 

p=0 


where 


>.(*,)= 


ext 


-p{j2rc$p } 

exp{j27i(j-l)dp} 


In  the  above  Equations,  the  following  quantities  are  defined  as 

GJp=^\)cos(eo)sin(<|)p -y)  ;  p=0,l,---,(Nc -l) 

^p=^-cos(e0)sin(<|)p)  ;  p=0,l,---,(Nc -l) 

4>p  =— 7C-t-pA<t) 
rc2+Hp(Hp+2rE)l 


0ft=-sin 


A<()= 


2rt(H,+rE) 
n 


(N,-l) 

4v. 


t)  = 


X,f, 


cx  PRF 


(45) 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


The  power  of  the  p-th  clutter  return  and  the  clutter  to  noise  ratio  (CNR)  are  defined  as 

(52) 


and 


.  /.  A  ^  PtTuA,2cCcGt(())p,eo)gr(<l)p>Oo) 
(47t)2N0Lsrc4 


(53) 


respectively. 


qc=s0WcLc 
so  =Yc  sin(\|/) 


H„r 


\|r=sm 


-i  i  **p 


1  +  - 


2r„ 


2rB 


(54) 

(55) 
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(56) 


(57) 


2nrc 

IT 


L  = 


Pr  _ 


— 7 — r=— — “7 — \  if  Prtan(\|0<rctan(pJ 

cos(\i/J  2fBcos(yc) 

—  tan^le $  if  pr  tan(\jfc  )>rc  tan(Ple ) 


[  sin(vc) 


Pr  = 


2f„ 


g.Mo^GoG.MoIs.Mo)- 

For  uniform  array  pattern,  Ga  (<j)p>  ©0) 1S  given  by 


o.MoW 


sin|^-[sin((t)p )-  sin((J)0  )]] 


J  sinj^-  [sin(<I)p )- sin((J>0  )]j 
sinjM  [sin(<t)p )-  sin  (<|>0  )]| 


for  -90°<((>p  <90° 


J  sin|^  [sin(<t>p )-  sin(<>0  )]1 


* 


for  90°<<J)p  <270° 


For  a  Dolph-Chebychev  array  pattern,  Ga  (<J>p  ,0O )  is  given  by 

-jj|FMfs0)®u]2  for  -90°<<{>p  <90° 

Gb^lFfeIfc)®u]2  for  90°<(})p  <270° ' 


Ga(4>P.e0)=| 


(58) 


(59) 

(60) 


(61) 


(62) 


where 

F  is  the  discrete  Fourier  transform  (DFT)  operator,  ®  is  the  Shur  (element  by  element) 
product,  and  u  is  the  J-dimensional  vector  of  Dolph-Chebychev  weights. 


?j  (f  s0  )— 


ex 


p{j27tfs0 } 


_exp{j2jt(J  -  l)fs0}J 


(63) 


with 
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(64) 


fso  =  (d/X)c°s(<]>o  )sin(<t)0 ) 


gc  (<t>p ,  0O )  is  the  element  power  pattern  given  by 

(  )  lC0S  2(<))p)  for  — 90°<<J)p  <90° 

g e  \t  p  ’  oj=|GbCos2(’(j)J  for  90°<(J>p  <270° 


gr(<l)p-eo)=ge(<l>p>eo)  is  the  receive  element  power  pattern. 


(65) 


4.3.5  Hybrid  Clutter  Model 

In  this  section,  the  hybrid  clutter  model  is  described.  This  model  is  referred  to  as  the 
hybrid  due  to  the  fact  that  it  combines  the  inherent  simplicity  of  the  representative  model 
with  the  sophistication  of  the  physical  model.  To  this  end,  let  the  (JxJ)  clutter  covariance 

matrix  Rc(m)=  E^(n)cH(n  -  m)}  be  given  by 

)exp{j2jtmS 

R»=°?  — -  -  ;  m>0  (66) 

EF.W 

p=0 

where  m  is  the  temporal  lag  and  the  vector  e,  (dp )  is  defined  as 

e,M= 


ext 


p{j27tT>p  } 
exp{j2rc(j  -  l)dp } 


(67) 
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5 

Estimation  Methods 


In  this  section,  we  describe  the  estimation  methods  used.  To  achieve  this,  we  first 
show  how  we  utilize  the  analytic  expressions  developed  earlier  in  a  process  synthesis 
procedure.  This  will  involve  generation  of  both  clutter  and  partially  correlated  signal  and 
interference  data  as  observed  from  a  phased  array  radar.  Let  xk  m(n)  be  a  zero-mean 

Gaussian  complex-valued  sample  representing  clutter,  target,  or  an  interference  signal 
observed  from  the  m-th  channel  corresponding  to  the  k-th  range  cell  at  time  n.  Consider  the 
space-time  data  vector  at  range  cell  k. 


'xk(l)‘ 

x.  - 

where  each  entry  xk(n)  is  a  Jxl  spatial  vector  given  by 


(68) 


xk 


xk,i(n) 

xk,2(n) 


x 


k,J 


and  J  and  N  are  the  number  of  channels  and  pulses,  respectively. 


Consider  now  the  JNxJN  Hermitian  block  Toeplitz  correlation  matrix 


R* 


Rk(-*) 

Mo) 


R,(N)  R,(N-1) 


R„(-N  +  l)' 

Rk(-N  +  2) 
Rk(0)  . 


(69) 


(70) 


where  each  JxJ  block  of  Rk  can  be  expressed  as 

R„(n)  R«(“)  —  R»(n) 

R21(n)  R22(n)  •••  R2J(n) 


Rk(n)  = 


[Ri.W  Rn(n)  -  R„(n)J 


;  n=0,l,...,(N-l). 


(71) 
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Next,  we  discuss  two  methods  for  data  generation. 


5.1  Block  Procedure  for  Clutter  Generation 

The  block  procedure  for  generating  clutter  and  partially  target  or  interference 
processes  is  very  useful  in  generating  independent  identically  distributed  (iid)  radar  data.  It 
is  also  one  of  the  fastest  data  generation  technique  for  quick  check  up  and  analysis.  Given  a 
known  clutter  covariance  matrix,  the  objective  is  to  factorize  this  matrix  using  techniques 
such  as  Cholesky  decomposition,  an  LDL  decomposition,  an  eigen-decomposition  or  a 
singular  value  (SVD)  decomposition.  Whereas  some  techniques  have  more  merit  than  others 
do,  we  have  shown  through  our  computer  simulations  that  the  SVBD  is  the  most  robust.  This 
in  turn  leads  to  more  stable  results.  Once  the  factorization  is  complete,  driving  a  white  noise 
process  through  this  system  generates  radar  data.  The  specific  steps  involved  in  this 
technique  are  described  below: 


1 .  Given  the  elements  of  the  JNxJN  clutter  correlation  matrix  at  range  k,  RCK ,  we  form  the 


SVD  of  Rck  , 


Rck  =  USUH  =  (US1/2)(US,/2)H . 
2.  Generate  a  data  vector  Xk  as 

xk=(us1/2)wk. 


(72) 

(73) 


where  Wk  is  a  JNxl  Gaussian  white  noise  vector  with  zero-mean  and  unit  variance.  In 
vector  form,  Equation  (49)  can  be  written  as 


Xk  = 

^3 

x\  XI 

1 _ 

=  (us,/2)wk 

=  (usI/2' 

w  w' 

...  ^  rr 

S  S 

i 

JU(N). 

_wk(N)J 

where 


xk,(n)l 


xu(n) 
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and 


wki(n) 

(  \  i  wk2(n) 
wk(n)  = 

wu(n)J 

The  elements  of  the  concatenated  JNxl  vector  Xk  consists  of  the  processes  containing  the 
correlation  specified  by  the  clutter  correlation  matrix  R* . 


5.2  Time  Series  Procedure  for  Clutter  Generation 

Using  this  technique,  the  clutter  signal  returns  are  generated  as  multi-channel  auto¬ 
regressive  (AR)  processes  [].  In  what  follows,  we  give  a  brief  description  of  the  proposed 
procedure. 

1 .  The  order  P  of  the  AR  process  is  selected  based  upon  a  specified  tolerance  for  fitting  the 
desired  spectrum. 

2.  The  elements  of  a  JNxJN  matrix  R  k  (this  matrix  is  a  subset  of  the  matrix  R  k )  are 
specified  using  the  analytic  expressions  developed  in  the  previous  sections. 

3.  Using  the  block  Toeplitz  form  of  Rk ,  the  multi-channel  Yule-Walker  equations  specified 
by  Equation  (51)  are  solved  using  the  Levinson-Durbin  recursion  to  determine  the  matrix 
coefficients  A”  (p)  and  the  Hermitian  driving  white  noise  covariance  matrix  [Zu  J1  for  p 
=  1,2,...,  P. 

A"Rl=ts.f  [0]  -  [0]}  (75) 


where 


AH=[l  AH(l)  Ah(2)  •••  ah(p)J 
[lJ=E|uuBJ=|xa 


'u,(n)' 

u2(n) 

_Uj(n)_ 


4.  Perform  the  SVD  of  the  matrix  [Eu  J1 ;  i.e., 


(76) 

(77) 

(78) 
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(79) 


[su]cH  =usuH  =(us,/2)(us1/2)h 

5.  These  matrices  are  now  used  in  the  generation  of  c(n)  using 

c(n)  =  -XAcH(p)c(n-p)+ uc(n)  (80) 

P=i 

where 

s»=  (u.sfkW  <81) 

and  vc(n)  is  a  white  noise  vector  whose  covariance  matrix  is  the  identity  matrix.. 
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6 

Examples 


In  this  section,  we  give  some  examples  showing  how  the  proposed  algorithms 
perform.  Note  also  that  part  of  this  work  set  the  stage  for  a  more  sophisticated  MATLAB- 
based  tool  referred  to  as  the  Air  Force  Research  Lab  (AFRL)  Multi-Channel  Signal 
Processing  System  (MCSPS).  This  tool  allows  the  researcher  to  have  a  complete  radar 
system  analysis,  starting  with  a  thorough  and  detailed  signal  generation,  choice  of  signal 
processing  technique  and  diagnostics.  It  is  the  first  part  of  the  algorithm  that  we  will  be 
using  in  this  section. 

To  simplify  the  analysis,  the  scenario  used  consists  of  J  =  14  channels  and  N  =  16 
pulses.  In  all  examples,  additive  white  Gaussian  noise  is  assumed  to  be  present.  Two  cases 
are  considered.  In  addition  to  the  clutter,  a  target  is  always  present,  but  no  jammer  will  be 
present  in  the  first  case  whereas  the  second  case  has  a  jammer  included  in  the  simulation.  In 
all  cases,  the  target  is  assumed  to  be  located  at  an  azimuth  angle  of  -30°  and  elevation  angle 
of  0°.  Moreover,  it  is  assumed  to  have  a  0.5  normalized  Doppler  and  a  Signal  to  Interference 
plus  Noise  ratio  (SINR)  of  10  dB.  In  all  cases  361  clutter  patches  are  assumed  to  make  up 
the  clutter  ring.  The  parameters  used  in  the  simulation  will  be  described  for  each  case.  We 
start  with  the  representative  model. 

6.1.  Representative  Model 

As  was  described  earlier,  this  model  generates  data  following  a  very  simple 
covariance  matrix  formulation.  In  the  following,  we  plot  the  angle/Doppler  map.  The 
scenario  consisted  of  clutter  and  a  single  target.  The  shaping  parameters  and  |is  have  been 
varied  so  as  to  generate  different  clutter  shapes. 
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Figure  9.  Angle/Doppler  map,  p,  =  0.9999,  p,  =  0.75. 
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Figure  10.  Angle/Doppler  map,  p*  =  0.75,  p  =  0.9999. 
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Figure  1 1 .  Angle/Doppler  map,  m  =  0.5,  p»  =  0.5. 


Figure  12.  Angle/Doppler  map,  p«  =  0.5,  p,  =  0.5, 1  jammer. 


Figure  9  shows  the  case  of  highly  correlated  clutter  (in  time),  whereas  the  spatial  correlation 
is  moderate.  This  leads  to  a  narrow  clutter  beam  in  the  angle/Doppler  map.  Figure  10  is 
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similar  to  figure  9.  The  values  of  (it  and  Li,  have  been  switched  and  therefore,  the  clutter 
ridge  shows  up  in  the  other  opposite  direction.  Since  the  target  remained  in  the  same 
location  (angle  and  Doppler),  whereas  it  was  in  the  clear  in  figure  9,  it  now  competes  with 
very  strong  sidelobe  clutter.  When  the  values  of  |it  and  (is  have  been  assumed  to  be  similar 
and  equal  to  0.5,  the  resulting  clutter  spectrum  resembles  that  of  a  bell  (see  figure  1 1). 

In  figure  12,  a  jammer  was  added  to  the  interference  and  shows  at  all  Doppler,  with  a 
normalized  spatial  frequency  of  0.2. 

6.2.  Physical  Model 

The  following  examples  are  taken  from  the  physical  model. 


-0.4  -0.2  0  0.2  0.4 

Figure  13.  Angle/Doppler  map,  v„  =  50. 
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Figure  15.  Angle/Doppler  map,  vp  =  50,  crab  angle  =  10°. 
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In  figure  13,  the  platform  velocity  was  assumed  to  be  50  m/s.  Along  with  other  parameters, 
this  leads  to  clutter  ridge  along  the  diagonal  in  the  angle/Doppler  map.  If  this  speed  is 
doubled  (100  m/s),  this  situation  creates  ambiguities  as  can  be  seen  in  figure  14.  In  both 
cases,  the  crab  angle  was  assumed  to  be  0°.  Figure  15  shows  the  case  when  the  crab  angle  is 
equal  to  10°;  this  leads  to  a  curved  clutter  ridge.  Again,  if  a  jammer  is  added  to  the  scenario, 
it  will  appear  all  Dopplers  at  a  specific  angle,  as  is  the  case  in  figure  16.  The  examples 
shown  above  are  simply  an  illustration  of  the  way  different  spectra  could  be  generated. 


6.2  Hybrid  Model 

In  this  section,  we  present  examples  taken  from  the  hybrid  model.  In  this  case  the 
clutter  ridge  slope  is  taken  to  be  and  a  temporal  clutter  one  lag  of  0.9  is  also  assumed.  Figures 
>  17-19  show  the  spectrum  of  the  described  scenario  as  a  function  of  |it. 
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Figure  18.  Angle/Doppller  map,  p,  =  0.9999. 
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Figure  19.  Angle/Doppler  map,  p,  =  0.5. 
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7 

Conclusions 


In  this  short  report,  we  described  three  techniques  for  data  generation.  These 
techniques  have  successfully  been  implemented  into  the  AFRL-MCSPS  simulation  and 
analysis  tool.  Engineers  at  AFRL  are  developing  and  studying  advanced  signal  processing 
techniques  such  as  STAP,  using  this  tool.  The  way  data  signals  are  generated  and  the 
confidence  in  the  approach  used  are  very  essential  for  a  good  and  thorough  performance 
evaluation.  For  these  reasons,  we  have  decided  to  include  all  three  techniques  as  part  of  the 
signal  generation  capability  of  the  tool.  In  section  6  of  the  report,  we  have  given  example  of 
spectra  generated  through  the  tool,  using  the  proposed  approaches.  We  have  shown  that 
these  spectra  are  representative  of  actual  airborne  data. 

Once  the  covariance  matrix  has  been  computed  using  any  of  the  three  proposed 
techniques,  space-time  data  is  generated  using  a  matrix  factorization  technique.  We  have 
tried  several  techniques  such  as  the  Cholesky  decomposition,  the  SVD  decomposition,  the 
LDL  decomposition,  the  eigen-decomposition,  and  the  square  root  of  the  matrix  (MATLAB 
function).  Of  all  the  techniques  mentioned  above,  we  found  that  the  SVD  performed  the  best 
in  terms  of  robustness  to  numerical  inaccuracies  as  well  as  speed  of  execution.  We  have  then 
decided  to  implement  only  this  technique  for  data  generation. 

Present  and  future  work  involves  developing  more  sophisticated  target  signal  and 
jammer  signals.  This  will  be  described  in  future  reports. 
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